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ABSTRACT 

We explore how the star formation efficiency in a protocluster clump is regulated 
by metallicity dependent stellar winds from the newly formed massive OB stars (M^ ^ 
5 Mq) on their main sequence. The model describes the co-evolution of the mass 
function of gravitationally bound cores and of the IMF in a protocluster clump. Dense 
cores are generated uniformly in time at different locations in the clump, and contract 
over lifetimes that are a few times their free fall times. The cores collapse to form 
stars that power strong stellar winds whose cumulative kinetic energy evacuates the 
gas from the clump and quenches further core and star formation. This sets the final 
star formation efhciency, SFEf. Models are run with various metallicities in the range 
Z/Zq — [0.1, 2]. We find that the SFEf decreases strongly with increasing metallicity. 
The ^FSj-metallicity relation is well described by a decaying exponential whose exact 
parameters depend weakly on the value of the core formation efficiency. We find that 
there is almost no dependence of the S'FiJy-metallicity relation on the clump mass. 
This is due to the fact that an increase (decrease) in the clump mass leads to an 
increase (decrease) in the feedback from OB stars which is opposed by an increase 
(decrease) in the gravitational potential of the clump. The clump mass-cluster mass 
relations we find for all of the different metallicity cases imply a negligible difference 
between the exponent of the mass function of the protocluster clumps and that of 
the young clusters mass function. By normalizing the SFEs to their value for the 
solar metallicity case, we compare our results to S'Fi?— metallicity relations derived 
on galactic scales and find a good agreement. As a by-product of this study, we also 
provide ready-to-use prescriptions for the power of stellar winds of main sequence OB 
stars in the mass range [5, 80] Mq in the metallicity range we have considered. 

Key words: galaxies: star clusters - Turbulence - ISM: clouds - open clusters and 
associations 



1 INTRODUCTION 

1.1 OBSERVATIONS OF THE SFE IN 

STELLAR CLUSTERS AND IN GALAXIES 

One of the most essential quantities that regulates the dy- 
namical evolution and chemical enrichment of stellar clus- 
ters, the interstellar medium, and galaxies is the star forma- 
tion efficiency (SFE) (e.g., Geyer & Burkert 2000; Boissier et 
al. 2001; Krumholz & McKee 2005; Dib et al. 2006; Parmen- 
tier & Fritze 2009; Dib et al. 2009). The SFE is commonly 
defined as the fraction of gas which is converted into stars 
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in a system of a given mass, be it a protocluster molecular 
clump, an entire giant molecular cloud (GMC), or a galaxy. 
In a non isolated system, such as a protocluster clump em- 
bedded in a filamentary structure from which the clump con- 
tinues to accrete, the final star formation efficiency, SFEf, 
is difficult to define. It can be approximated by: 

SFEf = SFEit = t..,) = . (t..,) 

where t^xp is the epoch at which the gas is expelled from 
the protocluster region, or otherwise depleted. M duster (texp) 
is the final stellar mass, Mgas,i is the initially avail- 
able reservoir of star forming gas, and Mgas,acc(texp) ~ 



2 S. Dib et al. 



j^texp M^i^^^^(t ) dt is the amount of gas the clump has ac- 
creted from the filament and from which also stars form, 
with Mci^acc being the clump's time dependent accretion 
rate. For a clump with low levels of accretion or for which 
star formation started after the bulk of the gas in the fila- 
ment have been accreted onto it, Eq.[T]is reduced to: 



SFEf 



Mr.. 



(2) 



where Mdump is the mass of the protocluser clump, and 
Mciuater is the mass of the stellar cluster. However, in the ob- 
servations, it is difficult to estimate what the original clump 
mass was, once a large fraction of the gas has been expelled 
from the protocluster region. Nevertheless, the SFE can be 
measured in embedded or semi-embedded clusters and is 
usually defined as being: 



SFEot 



Mri-u 



+ Mgas ' 



(3) 



where Mgas is the mass of the star forming gas in the proto- 
cluster region. Measured SFEs of nearby embedded clusters, 
using Eq.Oyield values that fall in the range 0.1-0.4 (Cohen 
& Kuhi 1979; Wilking & Lada 1983; Rengarajan 1984; Wolf 
et al. 1990; Pandey et al. 1990; Lada et al. 1991a,1991b; 
Warin et al. 1996; Olmi & Testi 2002). These SFE values 
are larger than those obtained for entire GMCs which are 
observed to fall in the range of a few percent (e.g., Duerr 
et al. 1982; Myers et al. 1986, Fukui & Mizuno 1991; Evans 
et al. 2009). On galactic scales, the star formation efficiency 
is usually defined as being the inverse of the molecular gas 
consumption time, and is given by: 



SFE„ai = 



SFRna 



gal 



M, 



H2 



(4) 



where SFRgai is the galactic star formation rate (in M0 
yr~^) and the mass of the molecular hydrogen gas. In 
some studies, is replaced by Mgaa which is the total 

mass of the gas in the galaxy, or by the mass of the H I 
gas (e.g., Schiminovitch et al. 2010). In large spiral galax- 
ies SFEgai ~ 0.5 X 10~^ yr~^ (e.g., Kennicutt 1989,1998, 
Murgia et al. 2002; Leroy et al. 2008). The local SFE in 
galaxies is found to depend on galactic radius (e.g., Zasov 
& Abramova 2006; Leroy et al. 2008). Among spiral galax- 
ies, the global SFEgai is found not to depend clearly on 
the Hubble type and to depend weakly on the galactic en- 
vironment (e.g.. Young et al. 1986; Rownd & Young 1999). 
However, Young et al. (1996) report that the SFEgai of Ir- 
regular galaxies (which are usually less metal-rich) is higher 
that that of normal spirals (Fig. 11 in their paper). Boissier 
et al. (2001) found no dependence of the SFEgai on some 
global galactic properties such as the circular velocity, the 
H band luminosity of the galaxies, and their B — H colour 
index. In Fig. [T] we have compiled the global SFEgai for 
a few selected nearby galaxies as a function of their global 
metallicity. The trend in Fig. [T] is strongly suggestive of an 
SFEg which increases with decreasing metallicity. Interest- 
ingly, a recent study by Mannucci et al. (2010) using SDSS 
data for nearby galaxies in combination with a sample of 
higher redshift galax;ies confirmed the existence, for a given 
galactic stellar mass, of a dependence of the SFRgai and 
of the specific star formation rate {SFRgai per unit stellar 
mass, SSFRgai) on metallicity. Their results show that the 
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Figure 1. Global galactic star formation efficiency SFEgai as a 
function of the galactic global metallicity. The data for M33 and 
NGC 6822 are from Gratier et al. (2010a,b), for the SMC from 
Leroy et al. (2006), and the average spiral value from the sample 
of Murgia et al. (2002). 



SFRgai and the SSFRgai increase with decreasing galactic 
metallicities. A similar result was also pointed out by Ellison 
et al. (2008) and is also visible when plotting the SSFR as 
a function of metallicity using the data of Lara-Lopez et al. 
(2010). 

1.2 REGULATION OF THE SFE IN STAR 
FORMING REGIONS 

Several physical processes regulate the core formation effi- 
ciency (CFE) within a star forming molecular clump/cloud, 
and as a consequence, the SFE. Supersonic turbulence 
that is ubiquitously observed in molecular clouds can, if 
driven, provide support against the global collapse of the 
clump/cloud, or at least delay it (e.g., Vazquez-Semadeni 
& Passot 1999). However, on scales smaller than the energy 
injection scale, but larger than the sonic scales in the cloud, 
supersonic turbulence produces local compressions (cores) 
(e.g., Padoan 1995; Vazquez-Semadeni et al. 2003) of which 
a fraction can be 'captured' by gravity and proceed to col- 
lapse into stars (e.g., Vazquez-Semadeni et al. 2005a; Dib 
et al. 2007a; Dib & Kim 2007; Dib et al. 2008a; Dib et al. 
2010a). Work by Clark & Bonnell (2004) and Klessen et 
al. (2000) suggested that the CFE decreases with increas- 
ing turbulent energy and decreasing energy injection scale. 
Krumholz & McKee (2005) formulated an analytical theory 
in which the star formation rate per free-fall time, SFRff, 
is shown to decrease with an increasing sonic Mach number 
and an increasing virial parameter (see also the simulation of 
Clark & Bonnell 2004). Padoan & Nordlund (2011) argued 
that the SFRff increases with increasing Mach number and 
that it also depends weakly on the ratio of the mean gas pres- 
sure of the cloud to the mean magnetic pressure. A similar 
trend for the dependence of the SFRff on the sonic Mach 
number was also reported by Rosas-Guevara et al. (2010), 
who additionally showed that the mass contained in collid- 
ing streams from which protocluster clumps form is also a 
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key parameter in regulating the fraction of the gas that can 
become gravitationally bound. 

Magnetic fields also play an important role in determin- 
ing the fraction of gravitationally bound gas in star form- 
ing clouds/clumps. Results from both ideal and non-ideal 
MHD simulations show that stronger magnetic fields (in 
terms of magnetic criticality) lower the rate of dense core 
formation in a star forming molecular clump/cloud (e.g., 
Vazquez- Semadeni et al. 2005b; Price & Bate 2008; Dib et 
al. 2008a; Li et al. 2010; Dib et al. 2010a). Dib et al. (2010a) 
showed that the CFE per unit of the free-fall time of the 
cloud, CFEff, are of the order of ~ 6 % and ~ 33 % for 
clouds with mass to-magnetic flux ratios of fiB = 2.2 and 
8.8, respectively (with /is being normalized by the critical 
mass-to-flux ratio for coUajise). 

Albeit magnetic fields and turbulence regulate the rate 
at which gas in a cloud/clump is turned into dense cores, 
and subsequently into stars, the SFE in a star forming re- 
gion will continue to increase over time until either the star 
forming gas is exhausted or otherwise expelled from the 
cloud/clump. The role of stellar feedback in setting the fi- 
nal value of the SFE by totally removing the gas has been 
investigated by several authors. The role of protostellar out- 
flows has been studied theoretically by Adams & Fatuzzo 
(1996), Matzner & McKee (2000), Matzner (2002,2007) and 
numerically by Nakamura & Li (2005,2007), Li & Nakamura 
(2006) and Li et al. (2010). Although protostellar outflows 
seem to play an important role in generating a self-sustained 
turbulence in a protocluster forming region, they do not in- 
ject enough energy that can lead to the gas removal from the 
region (e.g., Nakamura & Li 2007). Supernova explosions are 
an efficient way of removing gas from the protocluster region 
(e.g., Parmcnticr et al. 2008; Baumgardt et al. 2008). How- 
ever, such explosions occur after a few million years from the 
time massive stars have formed (i.e., the hydrogen burning 
phase of a 9, a 20, and a 40 Mq solar mctallicity stars are 
> 22,7, and 4 Myr, respectively; Meynet & Maeder 2000). 
Thus, supernova may play a role in the removal of the gas 
only in regions where other gas removal mechanisms would 
be inefficient due for example to a very low core and star 
formation efficiencies. 

Another form of stellar feedback is associated with O 
and B stars, in their main sequence phase, and eventually 
beyond. OB stars emit UV radiation which ionizes the sur- 
rounding gas and heats it to temperatures of ~ 7000—10^ K. 
This warm and ionized bubble provides the environment in 
which particles accelerated from the stellar surface by inter- 
action with some of the stellar radiation propagate outwards. 
In the first few hundred years, the wind freely expands at 
the wind velocity until it shocks with the surrounding ma- 
terial and/or with the wind of a neighbouring massive star. 
Shocked material is brought to temperatures of 10® — 10^ K. 
The hot gas expands supersonically in the warm medium. 
When the densities of the swept-up material reach a high 
value, a cooling instability operates behind the shock front 
leading to the formation of a thin dense shell. Models of en- 
ergy driven wind bubbles which take into account the inter- 
action of the expanding bubble with a surrounding confining 
medium were developed by Castor et al. (1975), Weaver et 
al. (1977), Shun (1980), Koo & McKee (1992a,b) and Canto 
et al. (2000), and recently included in global molecular cloud 
simulations by Vazquez-Semadeni et al. (2010). In contrast. 



models developed by Chevalier & Clegg (1985) and Stevens 
& Hartwell (2003) ignored the surrounding material and de- 
scribed the evolution of a free-floating wind. McKee et al. 
(1984) pointed out that the hot gas pressure in the bubble 
can be reduced by leakage of energy in the form of radiation. 
Harper-Clark & Murray (2009) argued that since the winds 
are likely to be clumpy (e.g., Lepine & Moffat 2008; Prinja 
& Massa 2010), the pressure in the bubble is reduced by the 
escape of the hot gas rather than the escape of radiation. 

For star clusters whose masses are < 10" Mq (corre- 
sponding to luminosities of lO^'' ionizing photon s~^ for 
a fully sampled IMF with a Salpeter slope) , the direct con- 
tribution of radiation pressure in the evacuation of gas from 
around the star/star cluster is minimal (e.g., Mathews 1969; 
Gail & Sedlmayr 1979; Arthur et al. 2004; Henney 2007; 
Krumholz & Matzner 2009). For more massive clusters, like 
those found in the Antennae, M82, and potentially some 
of the Milky Way starburst clusters such as Arches, several 
authors have recently pointed out that direct radiation pres- 
sure might be the dominant mode of gas expulsion from the 
protocluster region (Krumholz & Matzner 2009; Fall et al. 
2010; Murray et al. 2010). 

1.3 AIM OF THIS WORK 

An entirely unexplored area in star formation theories is the 
dependence of the star formation efficiency in a protoclus- 
ter environment on metallicity. As stellar winds are strongly 
metallicity dependent (e.g., Vink et al. 2001; Bresolin & Ku- 
dritzki 2004), feedback by winds is expected to lead to vari- 
ations in the SFE in protocluster regions which are identical 
except in the metallicity of the gas. To investigate this ef- 
fect, we use a modified version of the model developed by 
Dib et al. (2010b) which describes the co-evolution of the 
mass function of gravitationally bound cores and the IMF in 
a protocluster region. In this model, gas is evacuated from 
the protocluster clump by the stellar winds of the newly 
formed massive stars. In order to focus on the effects of the 
metallicity dependent feedback from stellar winds, we ne- 
glect in this paper the evolution of the core mass function 
(CMF) under the cfl'ects of clump gas accretion by the cores 
and core coalescence. The effects of these processes have 
been explored by Shadmehri 2004, Basu & Jones 2004; Dib 
2007, Dib et al. 2007b, 2008b, 2010b; Veltchev et al. 2010. 
In the Dib et al. (2010b) model, dense cores form in the 
protocluster clump uniformly in time following a specified 
core formation efficiency per unit free-fall time of the clump, 
CFEff. The dense cores have lifetimes of a few times their 
free-fall times, after which they collapse to form stars and 
are removed from the CMF. Stellar winds from the newly 
formed massive stars (M* > 5 Mq) inject energy into the 
protocluster clump. In order to calculate the stellar wind 
luminosities, we use a state of the art modified version of 
the steUar evolution code CESAM (Morel 1997, Morel & 
Lebreton 2008, Piau et al. 2011) to calculate main sequence 
models for stars in the mass range [5-80] Mq and derive 
their basic stellar properties; i.e., luminosity, effective tem- 
perature, and radius. The stellar grids are calculated for a 
range of metallicitics between Z = 1/10 Zq and Z = 2 Zq. 
In a second step, we use the stellar atmosphere model of 
Vink et al. (2001) in order to calculate the mass loss rates 
and terminal velocities of the winds. 
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2 THE CLUMP AND CORE MODELS 

2.1 PROTOCLUSTER CLUMPS: 

OBSERVATIONS AND MODELS 

; Several studies have established that star clusters form in 
dense (> 10'^ cm~^) clumps embedded in a lower density 
parental molecular cloud (e.g., Lada & Lada 2003; Lada et 
al. 2010; Csengeri et al. 2011; Parmentier 2011). Saito et al. 
(2007) recently studied, using the C^*0 molecular emission 
line, a large sample of cluster forming clumps whose masses 
and radii vary between [15-1500] MQand [0.14-0.61] pc, re- 
spectively. Based on the sample of clumps of Saito et al. 
(2007), the mass-size, and velocity dispersion-size relations 

' were fitted by Dib et al. (2010b) and are given by: 

M CUT \ in3-62±0.14 r>2.54±0.25/ \ ir\ 

Mciump(MQ) = 10 Jic (pc), (5) 
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Figure 2. Time evolution of the effective temperature (top), lu- 
minosity (middle), and radius (bottom) for solar metallicity stars 
of a few selected masses evolving on the main sequence. The cal- 
culation have been performed using the stellar evolution code 
CESAM (Morel 1997; Morel & Lebreton 2008). The three quan- 
tities are shown in real values (right column) and in percentage 
of variation with respect to the initial value (left column). 



In §.[21 we describe the protocluster clump model, the 
dense core model, and the distributions of cores that form 
in the clump (i.e., the initial CMF). In §.|3] we describe the 
metallicity dependent stellar wind models and the deriva- 
tion of the stellar wind luminosities. The co-evolution of the 
prestellar core mass function and of the IMF is presented in 
detail for a fiducial case with solar metallicity in i;. 14.11 We 
also show in this section how the feedback from stellar winds 
helps define the final SFE in the protocluster region, SFEf. 
The dependence of the of SFE on metallicity is presented in 
ij. 14.21 for the fiducial clump mass. In ^. 14.31 we investigate 
the dependence of SF^z-metallicity relation on the clump 
mass. In §.[S1 we discuss the implication of our results on the 
difference between the exponents of the clump mass function 
and the young clusters mass function and in §.[6] we discuss 
the implications of our results for the SFE-metallicty rela- 
tion on galactic scales as well as their implications for the 
SFE in globular clusters. In §.[71 we summarize our results. 



and 

n„ -1\ , n0'45±0.08 r>0.44±0.14 / 

Vc(km s ) = 10 itc (pc). 



(6) 



In this work, we adopt a protocluster clump model that 
follows an density profile: 



1 Pc{r) = 



PcO 



1 + {r/Rcor ' 



(7) 



where Rco is the clump's core radius (core radius here stands 
for the central region of the clump), pco is the density at the 
centre. In the study of Mueller et al. (2002), the average 
slope of the density profiles of 51 star forming clumps is 
— 1.8 ± 0.4 (here we adopt a value of —2). For a given mass 
Mciump of the clump, the central density pcO is given by: 



PcO = 



M, 



clump 



(8) 



- arctan{Rc/Rco)] ' 

where Rc is the size of the clump. The temperatures of the 
cluster forming clumps are observed to vary between 15 and 
70 K (e.g., Saito et al. 2007). In order to further constrain the 
models and minimize the number of parameters, we relate 
the sizes of the protocluster clumps to their masses using the 
mass-size relation of Saito et al. In the absence of detailed 
information about the velocity dispersion inside the clumps 
in the Saito et al. (2007) study, we assume that the clump- 
clump velocity dispersion they derived (i.e., Eq [6} is also 
valid on the scale of the clumps themselves and of their 
substructure. 



2.2 THE PRESTELLAR CORE MODEL 

Whitworth & Ward- Thompson (2001) applied a family of 
Plummer sphere-like models to the contracting prestellar 
dense core L1554, which is representative of the population 
of gravitationally bound cores in clumps that are considered 
in this work. They found a good agreement with the obser- 
vations of L1554 if the density profile of the core has the 
following form: 



Pp(r-p) 



PpO 



[1 + {rp/Rpor? ' 



(9) 



where ppo and Rpo are the central density and core radius 
of the core, respectively. Note that the radius of the core, 
Rp, depends both on its mass and on its position within the 
clump. The dependence of Rp on r requires that the density 
at the edges of the core equals the ambient clump density. 
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i.e., pp{Rp) ~ pc{r). This would result in smaller radii for 
cores of a given mass when they are located in their inner 
parts of the clump. The density contrast between the centre 
of the core and its edge is given by: 



C(r) 



PpO 

Pc{r) 



PpO 
PcO 



1 + 



r 

RcO 



(10) 



Depending on its position r in the clump, the radius 

of the core of mass M, Rp, can be calculated as being 
Rp{r,M) = a(r) Rpo{r,M), where: 



Rpoir,M) = 



M 
27rpj,o 



1/3 



(arctan[a(r-)]-3-^) ''',(11) 



and with a(r) = (C(r)^''^ — 1)^''^. With our set of parame- 
ters, the quantity C^^^ — 1 is always guaranteed to be pos- 
itive. The value Rp{r,M) can be considered as being the 
radius of the core at the moment of its formation. The 
radius of the core will decrease as time advances due to 
gravitational contraction. We assume that the cores con- 
tract on a timescalc, tcont,p which we take to be a few 
times their free fall timescale tff, and which is parametrized 
by tcontAr,M) = v tff{r,M) = (37r/32 Gp-p{r, M))^/\ 
where G is the gravitational constant, u is a, constant ^ 1 
and pp is the radially averaged density of the core of mass 
AI, located at position r in the clump. The time evolu- 
tion of the radius of a core of mass M, located at posi- 
tion r in the cloud is given by a simple contraction law 
Rp{r,M,t) = Rp(r,M,Q) e~(*/*-°"'.p). Theoretical consider- 
ations suggest that tcont,p can vary between tff {v — 1) and 
tff (v = 10) with the latter value being the characteritic 
timescale of ambipolar diffusion (McKee 1989; Fiedler & 
Mouschovias 1992; Ciolek & Basu 2001). Both observational 
(Lee & Myers 1999; Jessop & Ward- Thompson 2000; Kirk ct 
al. 2005; Hatchell et al. 2007; Ward- Thompson et al. 2007) 
and numerical (Vazquez-Semadcni et al. 2005a; Galvan- 
Madrid et al. 2007; Dib et al. 2008c) estimates of gravi- 
tationally bound cores lifetimes tend to show that they are 
of the order of a few times their free-fall time, albeit decreas- 
ing (but still larger than one free-fall time) when cores are 
defined with increasingly higher density tracers/thresholds 
One important issue is the choice of the cores central 
density, ppo. In this work, we first assume that the minimum 
density contrast that exists between the centre of the core 
and its edge is of the order of the critical Bonnor-Ebert value 
and that is > 15. Secondly, we assume that the density con- 
trast between the centre and the edge of the cores depends 
on their masses following a relation of the type ppo oc M**. 
Thus, the density contrast between the centre and the edge 
for a core with the minimum mass we are considering, Mmin 
(typically Mmin = 0.1 M©) is 15, whereas for a more mas- 
sive core of mass M, the density contrast will be equal to 
15 X {M/MminY- Observations show that fx varies in the 
range [0 — 0.6] (e.g., Caselli & Myers 1995; Johnstone & 
Bally 2006). 



2.3 THE INITIAL CORE MASS FUNCTION 

We assume that dense cores form in the clump as a result 
of its gravo-turbulent fragmentation. We describe the mass 
distribution of cores formed at each epoch using the formu- 



lation of Padoan & Nordlund (2002). Thus, the local mass 
distributions of cores, N{r,M), are given by: 



N{r,M) dlog M = /o(r) M 



-3/(4-/9) 



Jo 



P{Mj)dMj 



dlog M, 



(12) 



where /3 is the exponent of the kinetic energy power spec- 
trum, Ek oc , and is related to the exponent a of the size- 
velocity dispersion relation in the clump with l3 = 2a + 1. 
The local normalization coefficient fo{r) is obtained by re- 
quiring that / j^^^"*"" N{r, M) dM = 1 in a shell of width dr, 
located at distance r from the clump's centre. P{Mj) is the 
local distribution of Jeans masses given by: 



P{Mj) dMj = 



2 Mfo , 

■™ = M,/exp 



1 /In Mj ■ 

2 I ^ 



dMj, 



where Mjo is the Jeans mass at the mean local density, 
and aa is the standard deviation of the density distribution 
which is a function of the local thermal rms Mach number. 
Therefore, the local distribution of cores generated in the 
clump, at an epoch r, N{r, M, r), is obtained by multiplying 
the local normalized function N{r, M) by the local rate of 
fragmentation such that: 



N{r, M,T)dt - 



CFEff{r)p,{r) 



dt 



< M > (r) tcont,p{r, M) tfs,c 



N{r,M), (14) 



where dt is the time interval between two consecutive 
epochs, < M > is the average core mass in the local distribu- 
tion and is calculated by < M >= J^"*"" M N{r, M, 0) dM, 
and CFEff is a parameter smaller than unity which de- 
scribes the local mass fraction of gas that is transformed 
into cores per free fall time of the protocluster clump, f//,ci. 
In the present study, we assume that CFEff is independent 
of r. 



3 FEEDBACK MODEL: METALLICITY 
DEPENDENT STELLAR WINDS OF OB 
STARS 

We assume that the formation of cores in the protoclus- 
ter clump, and consequently star formation, is terminated 
whenever the fraction of the wind energy stored into motions 
that oppose gravity exceeds the gravitational energy of the 
clump. Whenever this occurs (at epoch t = texp), the gas is 
expelled from the protocluster clump and star formation is 
quenched. Thus, at any epoch t < texp, gas is removed from 
the clump only to be turned into stars. We take into account 
the feedback generated by the stellar winds of massive stars 
(M* > 5 Mq). In order to calculate reliable estimates of the 
feedback generated by metallicity dependent stellar winds, 
wo proceed in two steps. 

In the first step, we use a modified version of the stellar 
evolution code CESAM (see appendix 1 in Piau et al. 2011) 
to calculate a grid of main sequence stellar models for stars 
in the mass range [5-80] Mq (with steps of 5 Mq) at vari- 
ous mctalficitics which arc Z/Zq = [1/10, 1/6, 1/3, 1/2, 1, 2]. 
Provided the solar photosphere metal repartition recently 
determined by Asplund et al. (2005) is considered, the heavy 
elements (i.e., elements heavier than hydrogen and helium) 



6 S. Dib et al. 




20 40 60 80 

Mass (Mg^i) 



Figure 3. Time averaged, over the first 5 X 10^ yr of their life on 
the main sequence, effective temperature (top), luminosity (mid- 
dle), and radius (bottom) for stars of masses between 5 and 80 
Mq and metallicities ranging from Z = 1/10 Zq to Z = 2 Zq. 

fraction results from our latest solar calibrated models (Piau 
et al. 2011) and yields Zq = 0.0138. In addition to the solar 
metallicity case, the sample of explored metallicities is se- 
lected such that it covers a broad range of observed metallic- 
ities in nearby galaxies. This range covers the characteristic 
metallicities of nearby galaxies such as the Large Magellanic 
Cloud and M33 (Z/Zq ~ 1/3 - 1/2 (e.g, Russell & Dopita 
1992; Braine et al. 2001; RoUeston et al. 2002; Smith et al. 
2002), the metallicity of the SmaU Magellanic Cloud (SMC) 
ZsMC = 1/8-1/5 Zq (Russel & Dopita 1992; Rolleston 
et al. 2003; Lee et al. 2005; Perez-Montero & Di'az 2005), 
and metallicities characteristic of metal rich environments 
such as those prevailing in starburst galaxies such as M 83 
or in the circumnuclear starburst regions in the Galaxy and 
in other galaxies, Z — 1.5 — 2 Zq (Zaritsky et al. 1994; Kob- 
ulnicky et al. 1999, Bresolin & Kennicutt 2002). We always 
consider the solar metal repartition of Asplund et al. (2005) 
as we deal with typical population I stars. 

The evolution of massive stars was followed using the 
CESAM code for IMyr (or slightly less in some cases), on the 
main sequence. The characteristic stellar properties, which 
are the effective temperature Te//, the luminosity L-^, and 
the stellar radius ii* were dumped every 5 x 10* yr. For ages 
< 1 Myr, all the stars, including the most massive ones, re- 
main on the main sequence. Fig. [2] displays the time evolu- 
tion of Te//, Z/*, and Ri, in the models for a few selected 



stellar masses for the solar metallicity case, both in physical 
units (right column), and in the percentage of variation with 
respect to the initial value (left column). These figures show 
that the time variations of these quantities are of the order 
of a few percent, increasing with time to about 10-20 percent 
at most for the most massive stars. Since we are interested 
in the early time evolution of proto-clusters in which winds 
are injecting energy, we average the stellar properties over 
the first 5 x lO'' yr and use the averaged stellar properties 
as characteristic values for stars on the main sequence. This 
procedure has been repeated for all other metallicities. The 
resulting characteristic stellar properties are shown in Fig. [3] 
for the various metallicity cases. One can notice that for a 
given mass over the considered mass range, the radius of a 
higher metallicity star is larger and its effective temperature 
smaller thus resulting in luminosity values that are almost 
independent of metallicity. 

In a second step, we use the grid of calculated time av- 
eraged stellar properties to evaluate, for the different metal- 
licity cases, the stellar mass loss rates and the power of 
the stellar winds. To that purpose, we use the results of 
the stellar atmosphere models developed by Vink et al. 
(1999,2000,2001). These models allow us to evaluate the stel- 
lar mass loss rate M*, as a function of the stellar mass M*, 
effective temperature T^ff, the stellar luminosity L*, the 
metallicity Z, and the ratio of the velocity of the wind at 
infinity to the escape velocity, Voo/vasc, using the following 
formulations (Vink et al. 2001): 

log A'L = - 6.697 (±0.061) 

+ 2.194 (±0.021) log(i*/10^) 

- 1.313 (±0.046) log(Af*/30) 

- 1.226 (±0.037) log (^ '''"^T" ^ 
+ 0.933 (±0.064) log(reff/40000) 

- 10.92 (±0.90) {log(reff/40000)}^ 
± 0.85 (±0.10) \og{Z/ZQ) 

for 27 500 < Teff 50000K, (15) 

and 

log AL = - 6.688 (±0.080) 

± 2.210 (±0.031) log(L*/10'^) 

- 1.339 (±0.068) log(A/*/30) 

- 1.601 (±0.055) log (^^^^y^) 

+ 1.07 (±0.10) log(rcff/20000) 
+ 0.85 (±0.10) \og{Z/ZQ) 

for 12500 Tcff sC 22500K (16) 

where L* and M, are in solar units and T^g is in Kelvin. In 
this range the Galactic ratio of Uoo/fosc ~ 2.6 (Vink et al. 
2001). If the values for Voa at other metallicities are differ- 
ent from these Galactic values, then the mass-loss rates can 
easily be scaled accordingly. The escape velocity of the gas 
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at the surface of the star can be easily calculated using the 
standard formula: 



2 G A'h 



eff 



1/2 



(17) 



where G is the gravitational constant and Mt^^eff is the ef- 
fective stellar mass Mt,^i.ff — M«(l — H), where H is the 
factor that takes into account the correction of newtonian 
gravity by the radiation pressure due to electron scattering 
and which is given by (e.g., Lamers & Leitherer 1993): 



H = 7.66 X 10""cr, — ^ 



(18) 



where CTe is the electron scattering coefficient, which is given, 
for compositions in which the number of Helium and Hy- 
drogen atoms follows n{He)/{n{He) + n{H)) ~ 0.1, by 
(7e = 0.34 cm^ g-i if Te// ^ 35000 K, = 0.32 cm^ g-i 
when 30000 K ^ T^// < 35000 K, and = 0.31 cm^ g~^ for 
Teff < 30000 K (Pauldrach et al. 1990). Vink et al. (2001) 
did not derive the values of v^o , therefore, we use instead the 
relations obtained by Leitherer et al. (1992) in which Voo is 
given by: 



log 



3 (kms ^ ) = 



1.23 



0.30 log ( -1 



+ (Me 
+ 0.64 log(re//) 



+ 0.13 log j (19) 

Fig. |4]displays the mass loss rates calculated for the var- 
ious metallicities, using equations 1151 191 for OB stars in the 
mass range [5, 80] Mq (triangles). In order to allow the calcu- 
lation of mass loss rates for a star of any given mass, the data 
points in Fig. |3]were fitted using fourth order polynomials of 
the form log(M*) = X^^^q 4 Ai Ml for each of the considered 
metallicities. The values of the coefficients Ai, are reported 
in Tab.[T] The power of the stellar winds is given by M^v^. 
This quantity is displayed in Fig. [5] for the models with dif- 
ferent metallicities and the data points are fitted with a 
fourth order polynomial log(M* vl^) = J2i^o^Bi Ml. The 
coefficients Bi for the models with different metallicities are 
summarized in Tab. [21 The M*t;^ — M* relations displayed 
in Fig. [S] allow for the calculation of the total wind energy 
deposited by stellar winds. The total energy from the winds 
is given by: 



E,, 



120 Mq / N{M,)M4M^)vl 



dM^, ] dt'. (20) 



I Ml,=5 Mq 

We assume that only a fraction of E^ind will be trans- 
formed into systemic motions that will oppose gravity and 
participate in the evacuation of the bulk of the gas from 
the proto-cluster clump. The rest of the energy is assumed 
to be dissipated in wind-wind collisions or escape the wind 
bubble. The effective kinetic wind energy is thus given by: 



E. 



k,wind 



K E. 



wind 1 



(21) 



where k is a quantity ^ 1. It is currently difficult to esti- 
mate K as its exact value will vary from system to system 
depending on the number of massive stars, their locations, 
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Figure 4. Stellar mass loss rates for stars in the mass range 5- 
80 Mq on the main sequence, and for various metallicities. The 
stellar mass loss rates have been calculated using the stellar char- 
acteristics (effective temperature, stellar luminosity and radius) 
computed using the stellar evolution code CESAM coupled to the 
stellar atmosphere model of Vink et al. (2001). Over-plotted to 
the data are fourth order polynomial fits. The parameters of the 
fit functions are summarized in Tab.[T] 



and their wind interactions. As a conservative guess for the 
fiducial model, we take k — 0.1. Ek,wind is compared at every 
timestep to the absolute value of the gravitational energy. 



Egrav, which is calculated as being: 



pc{r) r dr, 



(22) 



3 JO 
where pc is given by Eq. [7] 

4 RESULTS 

4.1 THE CO-EVOLUTION OF THE CORE 

MASS FUNCTION AND OF THE IMF: A 
FIDUCIAL MODEL 

Whenever a population of cores of mass M, located at a 
distance r from the centre of the clump has evolved for a time 
that is equal to its contraction timescale, it is collapsed into 
stars. Thus, the local number of cores of a given mass, at a 
given epoch r, is the sum of all the local populations of cores 
of the same mass that have formed at all epochs anterior or 
equal to the considered epoch with the additional step of 
subtracting from that sum the cores of the same mass that 
have readily collapsed into stars. The local populations of 
cores of various ages are evolved separately as they are each 
in a different phase of their contraction (and if accretion onto 
the cores were included as in Dib et al. 2010b, they will have 
different accretion histories and rates), and will collapse and 
form stars at various epochs. Thus, the total local number 
of cores of a given mass M, at a time t, is given by: 



Ar(r,M,f) = ;^iV(r,M,r„i). 



(23) 
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Figure 6. Time evolution of tiie pre-stellar core mass function (left), and stellar mass function (right) in the protocluster clump with 
the fiducial model parameters. The last time-step shown is 1.79 t//,c which corresponds to the epoch at which gas is expelled from the 
protocluster clump. 



We assume that only a fraction of the mass of a core 
ends up locked in the star. We account for this mass loss in 
a purely phenonienological way by assuming that the mass 
of a star which is formed out of a core of mass M is given by 
M* = ^M, where $ < 1. Matzner & McKee showed that ^ 
can vary in the range 0.25—0.75. In the absence of strong ob- 



servational or theoretical constraints, we shall assume that 
^ is independent of the core mass. Finally, we account for 
the possible modification of the stellar mass function by the 
effect of stellar winds at the high mass end. The variations 
in the IMF at the high mass end will be given by: 
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Figure 5. The power of the stellar winds, or wind luminosities, 
for stars in the mass range 5-80 Mq on the main sequence, and 
for various metallicities. The stellar mass loss rates have been 
calculated using the stellar characteristics (effective temperature, 
stellar luminosity and radius) computed using the stellar evolu- 
tion code CESAM coupled to the stellar atmosphere model of 
Vink et al. (2001). The values of v^a have been calculated using 
the derivation by Leitherer et al. (1992). Over-plotted to the data 
are fourth order polynomials. The parameters of the fit functions 
are summarized in Tab. |2] 
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Figure 7. Time evolution of the ratio of the effective kinetic 
energy generated by stellar winds and the gravitational energy of 
the protocluster clump. Time is shown in units of the protocluster 
clump free fall timescale t//,c- The horizontal dashed corresponds 



to Ef, . 



1 with K = 0.1. The full line correspond to 



the fiducial model with the CFEjj = 0.2 and the dashed and 
triple dot-dashed to cases with CFEff = 0.1 and CFEff = 0.3, 
respectively. Diamonds correspond to the epochs at which the 
gas is evacuated from the cluster in the three models and the 
processes of core and star formation are terminated. 
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Figure 8. Time evolution of the SFE in the protocluster clump. 
The full line corresponds to the fiducial model with the CFE ff = 
0.2 and the dashed and triple dot-dashed to cases with CFEff = 
0.1 and CFEff = 0.3, respectively. Time is shown in units of the 
protocluster clump free fall timescale iff ,,. Diamonds correspond 
to the epochs at which the gas is evacuated from the cluster and 
indicates the final value of the SFE, SFEf, in the three models. 



dN^{r,M,t) 



dt 



(r, M,t), 



(24) 



where Af* is the stellar mass loss rate given by Eqs. 11511161 
and A^*(r,M,t) is the local number of stars, at time t, of mass 

In this section, we describe in detail the results for a 
fiducial model. In this model, the mass of the clump is taken 
to be yiciump = 10^ Mq, the metallicity is Z = Zq, and 
The CFE is CFEff = 0.2, where tff^d is the free-faU time 
of the clump and is given by t//,c; ~ (37r/32 GpcY^^, and 
Pc the average density of the clump. The other quantities 
have been taken to be equal to the most commonly cited 
observational determinations and have been set, in the fidu- 
cial model as well as in all other models, to the following 
values: Rco ~ 0.2 pc, u — 3, fi = 0.2, k = 0.1. We also as- 
sume that ^ = 1/3 which implies that 1/3 of the mass of the 
dense core ends up in the final star mass. It is worth men- 
tioning here that for clump masses of the order of a few 10^ 
Mq or less and in which stars form with a global SFE of a 
few up to a few tens of percents, the role of direct radiation 
pressure will be negligible and stellar winds are likely to be 
the dominant mechanism for the expulsion of the gas from 
the protocluster region. In this the luminosities of 

OB stars are almost independent of metallicity (i.e., Fig.|3]), 
any role played by direct radiation pressure will be to accel- 
erate the expulsion of the gas from the protocluster region 
without modifying the dependence of the stellar winds on 
metallicity. 

Fig. [5] displays the time evolution of the CMF (left col- 
umn) and of the IMF (right column) in the fiducial model. 
Each population of cores formed at a given epoch is the 
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sum of all populations of cores formed at different positions 
in the protocluster clump. The feature that appears at ~ 6 
Mq is a result of this summation. By t ~ 0.3 tff^ci, the first 
stars form. In this model, since /i — 0.2 > 0, the most 
massive stars form first, as the most massive cores tend 
to be, at any given position in the clump, more centrally 
peaked and thus have shorter lifetimes. As time advances, 
the IMF becomes fully populated. On the other hand, the 
CMF ceases to evolve (i.e., there is no accretion or coales- 
cence in this model as in Dib 2007, Dib et al. 2007b; Dib 
et al. 2008b; and Dib et al. 2010a) as the numbers of cores 
that are newly formed at each position in the protoclus- 
ter clump is balanced by an equal number of cores which 
collapses and forms stars. The final IMF of the cluster is 
established at ~ 1.79 f//,c, which corresponds to the last 
epoch shown in Fig. [6] At this epoch. Fig. [7| shows that the 
ratio Ek,mind/Egrav rcaches unity (full line), and as a conse- 
quence, gas is expelled from the protocluster region. Fig. [S] 
displays the time evolution of the SFE in this model (full 
line). At t = texp = 1.79 tff^ci, the final value of the SFE is 
SFEf ~ 9.05 X 10~^ 

In this fiducial model, we have adopted a CFE value 
of CFEff — 0.2. This is an intermediate value between the 
values of ~ 0.06 and ~ 0.33 measured by Dib et al. (2010a) 
in numerical simulations of molecular clouds with two differ- 
ent degree of magnetization (moderately and strongly super- 
critical clouds, respectively). At first glance, it may appear 
that increasing the CFE by a given factor will lead to an 
increase in SFEj by approximately the same factor. How- 
ever, for a fully sampled IMF, a larger CFE value implies 
that a larger number of OB stars will be formed and deposit 
larger amounts of feedback by stellar winds in the proto- 
cluster region. This in turn leads to a faster evacuation of 
the gas and to a limitation of the SFEf. Fig. [7] shows that 
in two other models similar to the fiducial case but with a 
different CFE, the evacuation of the gas occurs faster for the 
higher CFE case (i.e., case with CFEff = 0.3) and slower 
for the low CFE case (i.e., CFEff = 0.1). The time evolu- 
tion of the SFE in these two additional models is compared 
to the fiducial case in Fig. [S] The final values of the SFE 
at t = texp are SFEf ~ 6.59 x 10*^, ~ 9.05 x 10"^, and 
~ 0.11 for the cases with CFEff — 0.1, 0.2, and 0.3 respec- 
tively. A factor of 2 and 3 increase in the CFE leads to an 
increase of the SFEf only by factors of ~ 1.37 and ~ 1.64, 
respectively. This implies a strong regulation of the effect of 
a varying CFE by stellar feedback on the resulting SFEf in 
protocluster clumps. 

Another quantity that can modify the value of the 
SFEf is the mass-radius scaling relation of the protocluster 
clump. In this section and in the following ones, we have 
used the mass-radius relation of Saito et al. (2007) which 
is based on C^*0 (1-0) line observations of massive clumps. 
With this relation, the radius of a 10^ Mq clump is ~ 3.5 pc 
(Eq. [S]). Any other putative mass-radius relation that would 
for example result in the decrease of the radius of the clump 
Rc by a factor of 2 (for the same mass and clump core-radius 
values), would lead to the reduction of the gravitational en- 
ergy by a factor of ~ 4 (i.e., pco increase by a factor of 2, 
Eq. [8]) and an increase in the effective wind energy by a fac- 
tor of ~ 16 (enhancement in the number of OB stars by a 
factor of 2, and the feedback energy dependence is oc Mq^, 
where Mqb is the total mass of OB stars). In this case, the 



gas will be expelled earlier and the value of SFEf will be 
smaller. There are however some indirect indications that 
the Saito et al. (2007) relations provide a good description 
of the scaling relations that connect the masse, radii, and 
velocity dispersions of massive protocluster clumps. As al- 
ready discussed by Dib et al. (2010b), the bulk of the star 
formation in the clump occurs in the inner regions (between 
3-5 Rco) (see Fig. 8 in Dib et al. 2010b). With the adopted 
value of Rco = 0.2 pc, the bulk of the stellar mass will form 
in a region of up to ~ 1 pc radius. Cas expulsion is conducive 
to cluster expansion (e.g., Bastian & Goodwin 2006). Clus- 
ter expansion after gas removal will be visible in the form 
of a deviation from the initial density profile in the outer 
regions of the cluster and would increase the observed clus- 
ter radius by a factor of 2-2.5 (i.e., for a cluster mass of 
Mci — 7.5 X 10* M0 after 20 Myrs and by a similar factor 
for a 5 X lO'' Mq after 40 Myrs ). This would lead the ex- 
panding clusters formed from clumps that follow the Saito 
et al. (2007) scaling relations to have radii ~ 2 — 2.5 pc. This 
agrees very well with the observed mean value of the stellar 
clusters radii in M51 obtained by Scheepmaker et al. (2007). 

4.2 THE METALLICITY-STAR FORMATION 
EFFICIENCY RELATION 

We now turn to the effect of varying the metallicity in the 
protocluster clump. Fig. [5] displays the time evolution of 
the ratio Et.wind/ Egrav in models similar to the fiducial 
case but with metallicities varying between Z = 0.1 Zq and 
Z = 2 Zq. All six models have the fiducial value of the CFE, 
CFEff = 0.2. In models with lower metallicities, the power 
of the stellar winds in weaker (i.e.. Fig. [5]) and the evacuation 
of the gas occurs at later epochs as compared to the higher 
metallicity cases. In the model with Z = 2 Zq, the gas is 
expelled from the protocluster region at t ~ 1.4 t//,c! while 
in the model with Z — 0.1 Zq, the gas expulsion is delayed 
until t ~ 4.5 t//,c! . For a given CFE, longer timescales imply 
a larger SFEf. Fig. 1101 displays the time evolution of the 
SFE in all six models. The curve is similar for all six models 
as the structural and dynamical properties of the models 
are similar (only the metallicty is different). The different 
symbols on the curve in Fig. \TU\ mark the epochs at which 
gas is expelled from the protocluster clump for the different 
metallicity cases and indicate the corresponding values of 
the SFEf in each model. These values are plotted versus 
metallicity in Fig. 1111 A clear trend is observed in which the 
SFEf increases with decreasing metallicity (diamonds). We 
fit the SFEf — Z points with the following functional form: 

SFEf = Cz e'^^°''i^) , (25) 

with Cz = 0.091 ± 6.8 x 10"'' and tz = 0.91 ± 7.9 x 10"^ 
We have also repeated the calculations at the various metal- 
licities using the additional values of CFEff = 0.1 and 
CFEff = 0.3. The dependence of the SFE on metallic- 
ity in these cases are also shown in Fig. 1111 As in the solar 
metallicity case, Fig. llll shows that the SFEf vary typically 
by a factor of ~ 1.3 — 1.4 and ~ 1.6— 1.7 for variations of the 
CFE by factors of 2 and 3, respectively. The curves corre- 
sponding to the cases with CFEff = 0.1 and CFEff = 0.3 
were also fitted with the functional form given in Eq. 1251 
The values of the fit parameters are summarized in Tab. |4] 
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Figure 9. Time evolution of the ratio of the kinetic energy gen- 
erated by stellar winds and the gravitational energy of the proto- 
cluster clump. Time is shown in units of the protocluster clump 
free fall timescale t//,c- The horizontal dashed corresponds to 



1 with 



0.1. Diamonds correspond to the 



epochs at which the gas is evacuated from the cluster in the three 
models. The mass of the clump in these models is 10^ Mq. 
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Figure 10. Time evolution of the SFE in the protocluster clump. 
Time is shown in units of the protocluster clump free fall timescale 
tff^c- The symbols correspond to the epochs at which the gas is 
evacuated from the cluster and indicates the final value of the 
SFE , SFEf, for the models with different metallicities. 



4.3 THE DEPENDENCE OF THE SFE-Z 
RELATION ON THE MASS OF THE 
PROTOCLUSTER CLUMP 

In the previous section, we have explored the dependence 
of the SFE on metallicity for a fiducial clump mass of 10^ 
Mq. For a given value of the CFEff, If the mass of the 
clump decreases (increases), the number of stars includ- 
ing OB stars decreases (increases) and the feedback from 
stellar winds decreases (increases). Simultaneously however, 
when the clump mass decreases (increases), the gravita- 



Figure 11. S F E f-yietallicity relations for clumps with the fidu- 
cial mass of 10^ and for different core formation efficiencies 
of CFEff = 0.1, 0.2, and 0.3. Over-plotted to the data are fit 
functions (Eq. [25j whose parameters are listed in Tab. (4] 




Figure 12. Time evolution of the ratio of the effective kinetic 
energy generated by stellar winds and the gravitational energy of 
the protocluster clump. Time is shown in units of the protocluster 
clump free fall timescale t//,c- The horizontal dashed corresponds 



to Ek 



nd/Eg 



1 with K = 0.1. Diamonds correspond to the 



epochs at which the gas is evacuated from the cluster in the three 
models. The mass of the clump in these models is 5 X 10'* Mq. 



tional energy of the clump decreases (increases) and this 
allows for a shorter (longer) delay before gas expulsion from 
the protocluster clump starts hence for the buildup of a 
smaller (higher) SFE. In order to explore the effects on the 
SFEf — Z relation of these two opposing dependencies, we 
have performed additional models with clump masses in- 
creased and reduced by a factor of 2 from the fiducial clump 
mass value. Fig. [12] displays the time evolution of the ra- 
tio Ek,wind/ Egrav in the model with Md = 5 x 10'' M0. 
In comparison to cases with the fiducial mass (i.e.. Fig. 
the evacuation of the gas from the protocluster region oc- 
curs on shorter timescales. In contrast, Fig. 1131 shows that 
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Figure 13. Time evolution of the ratio of the effective kinetic 
energy generated by stellar winds and the gravitational energy of 
the protocluster clump. Time is shown in units of the protocluster 
clump free fall timescale t//,c- The horizontal dashed corresponds 
to Ej. uiind/ ^grav = 1 with K = 0.1. Diamonds correspond to the 
epochs at which the gas is evacuated from the cluster in the three 
models. The mass of the clump in these models is 2 X 10^ Mq. 
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Figure 14. 5Fi?j^-Metallicity relations for clumps of different 
masses with the fiducial core formation efficiency of CFEff = 
0.2. Over-plotted to the data are fit functions (Eq. I25II whose 
parameters are listed in Tab. \5\ 



the evacuation of the gas from a clump with a mass twice 
as large as the fiducial mass (i.e., 2 X 10^ Mq) occurs on 
timescales (in units of the corresponding free-fall timescale 
of the clump) larger than those with the fiducial mass for 
the different metallicity cases. 

Fig. [14] compares the SFEf — Z relation for the models 
with the three clump masses. The effect of reducing (in- 
creasing) the amount of feedback from winds by decreasing 
(increasing) the clump mass is compensated by the effect of 
the reduction (increase) of the gravitational potential. This 
results in final SFE values at the different metallicities that 
are nearly identical and independent of the clump mass. The 



fit parameters for the data displayed in Fig. [14] are summa- 
rized in Tab. [5] This negligible dependence of the SFEf — Z 
relation on clump mass is expected to remain valid for stellar 
clusters originating from clumps whose masses are smaller 
than the ones considered here, at least when averaged over a 
large number of clumps. On the other hand, the case of more 
massive clumps which would result in stellar clusters whose 
masses are larger than those considered here would require 
taking into account the effect of direct radiation pressure 
(e.g., Krumholz & Matzner 2009). This is left to a subse- 
quent work. Note however that very massive clusters are 
only found in large numbers in starburst galaxies such as 
the Antennae or M82. 

The range of SFEf values we obtained as a function of 
the adopted ranges in metallicity and CFEff is 0.05 — 0.35. 
A number of studies which have measured the fraction of 
bound stars and cluster survival following an instantaneous 
gas expulsion of the gas from the protocluster region indi- 
cate that a minimum value of SFEf ~ 0.33 is necessary for 
cluster survival (Lada et al. 1984; Geyer & Burkert 2001, 
Boily & Kroupa 2003a,b; Fellhauer & Kroupa 2005; Baum- 
gardt & Kroupa 2007; see summary in Fig. 1 of Parmen- 
tier & Gilmore 2007jj. It is important to note that these 
simulations did not explicitly model the dynamics of the 
gas. Instead, gas expulsion was performed by a reduction of 
the cluster potential which mimics a gas depletion process 
over a given timescale. For the range of SFEs found in our 
work, only models in which gas was removed on timescales 
of the order of 10 times the dynamical timescale of the clus- 
ter resulted in bound clusters. In these latter models, stars 
were assigned initial velocity dispersions such that the en- 
tire system (gas-|-stars) is in virial equilibrium. This most 
likely contradicts the results of many gravo-turbulent star 
formation simulations in which stars form in dense and qui- 
escent regions of the clouds/clumps and in which most of 
the turbulent motions have been dissipated (e.g., Klessen et 
al. 2000, Dib et al. 2007a). In the framework of the theory 
of gravo-turbulent fragmentation of a cloud/clump, which 
we envision to be responsible for core formation in the pro- 
tocluster clump considered here, stars will inherit the low 
levels of velocity dispersions from their gaseous progenitors. 
Geyer & Burkert (2001) performed a direct Nbody-|-SPH 
simulations in order to model the behaviour of stars and gas 
more consistently while requiring that stars form with zero 
velocity dispersion. They found that this configuration leads 
to the formation of compact clusters and to a severe reduc- 
tion of the SFE required to form bound clusters (50 percent 
of bound stars for a SFE of 0.05 and close to 70 percent for a 
SFE of 0.1; Fig. 7 in their paper). Further work is obviously 
needed to further clarify the dependence of the stellar clus- 
ter survivability on the SFE using models that implement 
more realistic gas expulsion schemes. 



^ One way of increasing the final SFE, is by increasing the 
CFEff. We have adopted realistic values of the CFEff in the 
range 0.1 - 0.3. Adopting CFEff = 0.5 will shift the SFEs ob- 
tained at the different metallicities by a factor of ~ 1.5 from their 
current values with the CFEff of 0.2 — 0.3. However, the SFEf 
for the solar metallicity case will still remain close to ~ 0.15 
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Figure 15. Clump mass-embedded stellar cluster mass relation 
in models of different metallicities and for the fiducial CFE value 
of CFEff = 0.2. 



5 IMPLICATIONS FOR THE EMBEDDED 
CLUSTER MASS FUNCTION 

An interesting application of the results presented above is 
to explore the relationship between the exponent of the pro- 
tocluster clumps mass function and that of the mass func- 
tion of embedded clusters. Furthermore, it is interesting to 
explore to which extent this relation depends on metallicity. 
Fig. [15] displays the stellar clusters masses as a function of 
their parental gas clump masses for the various metallicities 
considered in this work. All cases presented here are for the 
fiducial CFE value of CFEff = 0.2. The M^ump - Monster 
relations in Fig. [15] can be approximated by power laws of 
the form Mduster = Cm If the clump mass func- 

tion follows a power law of the form dN/d Ad clump oc M^i^mp 
and the young clusters mass function another power law 
given by dN/dMduster oc Mciusten then the exponents 7, 5, 
and 77 are bound by the relation 5 = 7 (77 + 1) — 1. Tab. [S] 



summarizes the values of 7 for the Mdump — Mduster re- 
lations for the various metallicity cases. The derived values 
of 7, which are close to unity, suggest that there should ex- 
ist only a small difference between the exponent of all the 
mass functions of protocluster clumps and that of the young 
stellar clusters. In fact, if 77 falls in the range —1.80 to —2 
as is observed for young stellar clusters in the Galaxy and 
in other spiral galaxies (e.g., Elmegreen & Efremov 1997; 
McKee & Wilhams 1997; Zhang & FaU 1999; Dowell et al. 
2008; Larsen 2009; Chandar et al. 2010), then the derived 
values of S fall in the range [-1.753, -1.808] (for 77 = -1.75) 
and [-2.004, -2.077] (for 77 = -2), for metallicities in the 
range of Z/Zq — [0.1,2], respectively. 

In a recent paper. Fall et al. (2010) suggested that a mo- 
mentum driven feedback, of any kind, is sufficient to relate 
the exponent of the young clusters mass function (77 ~ —2) 
and the value of 5 ~ —1.7 that is quoted by the latter au- 
thors as being the exponent of the protocluster clumps mass 
function. In reality, there is little evidence that the exponent 
of S ^ —1.7 relates to the mass functions of gravitationally 
bound structures. Dib et al. (2008b) showed that the mass 



function of gravitationally unbound, turbulence dominated, 
cores/clumps is well described by a power law whose ex- 
ponent is ~ —1.7. A similar result has been obtained by 
Hennebelle & Chabrier (2008). Dib et al. (2008a) also found 
that the mass function steepens with the increasing den- 
sity threshold and therefore increasing level of the gravi- 
tational boundedness of the selected cores/clumps, reach- 
ing values of > —2 for populations of cores/clumps which 
are entirely gravitationally bound (i.e., this happens when 
cores/clumps are detected at threshold densities of > 10^ 
cm~^, Dib et al. 2007a,, Dib & Kim 2007). Using a similar 
argument based on the existence of a threshold density for 
star formation to occur in the dense regions of lower den- 
sity molecular clouds/clumps (as argued for by Dib et al. 
2007a), Parmentier (2011) showed that the explains the dif- 
ference between the slope of the mass function of low density 
clouds/clumps (~ —1.7) and that of high density protoclus- 
ter forming clumps —2). 



6 IMPLICATIONS FOR GALAXIES AND 
GLOBULAR CLUSTERS 

There is, to date, no known observationally derived SFEf — 
Z relation for individual stellar clusters that can be directly 
compared to our models. It is however possible to establish 
such relation on galactic scales. Fig.[T]clearly shows that the 
galactic SFE, SFEgat, decreases with increasing metallicity 
of the gas. The negligible dependence we find of the SFEf 
on the clump mass suggests that this SFEf — Z relation 
remains valid on galactic scales. This is true as long as star 
formation in clusters is the main culprit in establishing this 
relation and that other processes such as gas infall onto the 
galactic disk and galactic outflows play only a minor role. 
One way of comparing our SFEf — Z relation(s) to the 
galactic SFEgai — Z relations is to calculate the relative 
SFE variations with respect to a characteristic value. For our 
stellar cluster models, we chose this value to be the SFEf 
at solar metallicity and for fiducial CFE value of CFEff — 
0.2. For galaxies, we normalize the SFEgai by the 'average 
spirals' value. Fig. 1161 displays the relative SFE, SFErei as 
a function of metallicity expressed in 12 -I- log[0/H]. The 
normalized relation corresponding to the fiducial case with 
CFEff — 0.2 is fitted with the functional form given by: 



SFEr 



SFEf 



SFEf{Z = Zq, CFEff = 0.2) 



■ Cz,rel e 



(26) 



yielding Cz,rei = 1.007±7.5 x lO'^* and Tz,rei = 0.915± 
7.0 X 10~^. The normalized relation corresponding to the 
cases with CFEff = 0.1 and 0.3 yield fit coefficients and 
exponents which differ by about ~ 0.3 dex and 0.03 from 
the fiducial case, respectively. The detailed values of all fit 
parameters are reported in Tab. [7] 

In this work, we have explored the dependence of the 
final SFE on metallicity in the metallicity range Z/Zq — 
[0.1, 2]. The power of stellar winds is expected to decrease 
further at metallicities lower than Z/Zq — 0.1. Thus one 
may expect that for halo like metallicities which are those of 
old globular cluster, the SFEf will be substantially higher. 
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Figure 16. Relative SFE as a function of metallicity. The rel- 
ative SFE for the models is taken to be SFEjj/SFEj{Z = 
Zq, CFEff^ = 0.2). The relative SFE for the observational data 
is SFEgai/SFEgai{a,ver.spir.). The full line, dashed line, and 
triple-dot dash lines corresponds to the model with CFEjj =0.2, 
0.1, and 0.3, respectively. 

The globular clusters in the Galcixy and in nearby galax- 
ies are observed to be more massive than 'standard' open 
clusters (e.g., Mandushev et al. 1991; Fan et al. 2010). One 
explanation of this is that globular clusters formed in an 
environment of high density and low angular momentum. 
They are also thought to have been more massive at birth 
but were affected by evaporation processes since the time of 
their formation. Our results which suggest, if extrapolated 
to halo like metallicities, that the SFEf in globular clusters 
parental clumps is very high further strengthens the idea 
that the initial masses of globular clusters at the moment 
of their formation must have been many times higher than 
their currently measured mass (e.g., Schaerer & Charbon- 
nel 2011). It should be noted however that many globular 
clusters show signature of multiple stellar populations (e.g., 
Bekki 2011 and refereces therein). The gaseous ejecta from 
the first generation of AGB stars can be converted into new 
stars and this will alter the original SFEf — Z relations. In a 
future work, we plan to investigate the SFEf in clumps with 
halo like metallicities, paying particular attention to the ef- 
fect of initial conditions (i.e., the shape and characteristic 
mass of the CMF). 



7 SUMMARY 

We use semi-analytical modelling to study the effects of gas 
dispersal from a protocluster clump by stellar winds of vari- 
ous metallicities in the range Z/Zq = [0.1, 2]. In this model, 
populations of gravitationally bound cores form uniformly in 
time at different positions from the protocluster clump cen- 
tre by turbulent fragmentation. The cores collapse to form 
stars after evolving over their lifetimes which are taken to 
be a few times their free fall times. The population of newly 
formed OB stars power strong stellar winds whose cumula- 
tive energy counters gravity and acts to disperse gas from 
the protocluster clump and hence quench further core and 



star formation. In order to estimate the power of the stellar 
winds, we use the stellar evolution code CESAM to calcu- 
late the basic stellar properties (i.e., effective temperature, 
stellar radius, and luminosity) of main sequence OB stars in 
the mass range [5 — 80] Mq. These calculations were per- 
formed for the metallicities of Z/Zq = 1/10, 1/6, 1/3, 1/2, 1 
and 2. The latter quantities were then used in the stellar 
atmosphere model of Vink et al. in order to derive the mass 
loss rates and power of the winds. Our main results are the 
following: 

a) For a given core formation efficiency (CFE), ex- 
pressed in units of the protocluster clump free-fall time 
CFEff, the final star formation efficiency, SFEf, which 
is set at the epoch at which the gas is evacuated from 
the protocluster clump, decreases with increasing metallic- 
ity. With the fiducial value of CFEff — 0.2, we find that 
SFEf = Cz e-'°s(^/^o)/^2, with Cz = 0.091 ± 6.8 x 10"'' 
and Tz = 0.91 ± 7.9 x 10"'^. This result is a direct conse- 
quence of the fact that lower metallicity stars have weaker 
winds than their higher metallicity counterparts. Weaker 
winds delay the evacuation of the gas from the protocluster 
clump and lead to higher values of the final star formation 
efficiency. With the fiducial value of the CFE, the SFEf 
varies between a few percent (~ 0.06 for Z = 2 Zq) to a few 
tens of percent (~ 0.27 for Z = 0.1 Zq). 

h) For a given clump mass, the SFEf, depends weakly 
on the CFE. Variations of the CFE by factors of 2 and 3 
lead to variations of the SFEf by factor of ~ 1.3 and ~ 1.6, 
respectively. 

c) For a given CFE, the value of the SFEf depends very 
weakly on the mass of the protocluster clump. This is due to 
the fact that any increase in the clump mass which results 
in an increase in the number of OB stars, and hence feed- 
back from stellar winds, is balanced by the increase in the 
gravitational potential of the clump which acts to delay gas 
evacuation. We find that varying the clump mass by a factor 
of 2 around the fiducial value of 10^ Mq in the metallicity 
range we have explored changes the SFEf only by factors 
of a few percent. This indicates that if the star formation 
efficiency in clusters is the main driver of the galactic star 
formation efficiency-metallicity relation, then our result of 
SFEf = Cf e-'°s(^/^o)/^z^ with Cz ~ 0.1 and rz ~ 0.91 
remains valid on galactic scales. We show that, when the 
cluster and galactic SFEs are normalized to a characteristic 
value (i.e., the corresponding SFE at solar metallicity), they 
are nearly identical (Fig. [16]). 

d) We find that the exponent, 7 of the clump mass- 
stellar cluster mass relation, (i.e., Mdump oc A^Jj^^^g^) is 
close to unity and increases only very weakly with decreasing 
metallicity. This implies that the exponent of the young star 
cluster mass function and that of the gravitationally bound, 
protocluster clumps mass function are nearly identical. 

e) As a by-product of this study, we provide fit functions 
of the power delivered by stellar winds from OB stars on the 
main sequence in the mass range [5, 80] Mq and for various 
metallicities in the range Z/Zq = [0.1,2]. These relations 
can be easily implemented in other analytical or numerical 
works studying the effects of feedback from stellar winds on 
the interstellar medium in galaxies of various metallicities. 
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Table 1. Coefficients of the fourth order polynomial fit to the Af*-logio (A/*^) relation {log{M*) = 5I]i=o 4 applied to the data in 

Fig- Si- The stellar mass, Mi,, is in solar mass units, Mq, and the stellar mass loss rate, A/* is in Mq yr^-*^. 
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Table 2. Coefficients of the fourth order polynomial fit to the Af*-logio(A^* i)^) relation (log(Af* f^) = J^i=o 4 ^» applied to the 
data in Fig. [5}. The stellar mass, Af«, is in solar mass units, Mq, and the power of the stellar winds (i.e., energy per unit time) in A/* 
is in J s"^. 
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Rp 


Radius of the core 




Free-fall timescale of the core 


icont ,p 


Contraction timescale of the core 


V 


Ratio of the contraction timescale of the core to its free-fall time 


Stellar and Feedback variables 


Meaning of the variables 


Af* 


masses of stars 




fraction of the mass of the core that ends up in the stars 


Af* 


Mass loss rates from stars 




Power of the stellar winds 




time integrated wind energy from OB stars stars 




Fraction of the wind energies that opposes the total gas-|-stars gravity in the clump 


^ k,wind 


Effective wind kinetic energy 



Table 3. Main variables of the model. From top to bottom, the first panel describes the protocluster clump variables, the second panel 
the dense cores variables, the third panel the stellar and stellar feedback variables. 
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CFEff 0.1 0.2 0.3 

Cz 6.62 X 10-2 ± 6.13 X IQ-"' 9.11 x 10"^ ± 6.89 x 10""' 1.09 x IQ-i ± 7.30 x IQ-" 
Tz 9.32 x 10-1 ± 1.01 X 10-2 9.15 X 10-1 ± 7.97 X 10-^ 9.01 x IQ-^ ± 6.90 x IQ-^ 

Table 4. Coefficients and exponents of the relation between the final SFE of the cluster SFEj and the metallicity of the gas for 
the case with the fiducial case (with CFEff = 0.2) and 2 other models with different CFE values {CFEff = 0.1 and 0.3) i.e., 
SFEf = Cz e-l°s(2/^o)/^z (Eq.[25t. 



Mc,„„p 5 X 10-* M0 10^ M0 2 X 10^ M© 

Cz 8.73 X 10-2 ±4.40 X 10""' 9.11 x 10"^ ± 6.81 x 10-* 9.35 x IQ-^ ± 9.32 x 10"* 
TZ 8.82 X 10-1 ±4.95 x 10-3 9.15 x lO-i ± 7.97 x 10-^ 9.36 x lO-i ± 1.10 x lO'^ 

Table 5. Coefficients and exponents of the relation between the final SFE of the cluster SFEf and the metallicity of the gas for 
SFEf = Cz e-^°siZ/ZQ)/Tz (Eq.[25ll, for various clump masses and for the fiducial CFE of CFEff = 0.2. 
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Z = 1/10Zq Z = 1/6Zq Z = 1/3Zq Z = 1/2Zq Z = Zq Z = 2 Zq 

7 1.004 1.011 1.024 1.033 1.047 1.077 

Cm 0.259 0.186 0.114 0.085 0.052 0.027 

Table 6. CoefScients and exponents of the relation between the clump masses and corresponding stellar cluster masses in our models, 
^Icluster = Cm ^^clump cases with different metallicities. 
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Tz'rei 9.32 X 10~i ± 1.01 X 10^2 9 15 ^ iQ-i ± 7.97 x lO^^ g gi X 10-1 ± 6.90 x 10"^ 

Table 7. Coefficients and exponents of the relation between the final SFE of the cluster SFEf, normalized to the solar metallicity value 
for the fiducial case with CFEff = 0.2, and the metallicity of the gas for the case with the fiducial case (with CFEff = 0.2) and 2 other 
models with diflferent CFE values {CFEff = 0.1 and 0.3) i.e., SFEr^i = SFEf/SFEf{Z = Zq, CFEff = 0.2) = Cz e" i°s(^/^o)/^z 
(Eq.[26j. 
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